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Right-handed neutrinos offer us the possibihty of accommodating neutrino masses. In a 
supersymmetric model, this imphes the existence of right-handed sneutrinos. Right-handed 
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■ sneutrinos are expected to be as hght as other supersymmetric particles if the neutrinos are 

! Dirac fermions or if the lepton-number breaking scale is at (or below) the supersymmetry (SUSY) 

I ' 

1' breaking scale, assumed to be around the electroweak scale. Depending on the mechanism of SUSY 

• breaking, the lightest right-handed sneutrino may be the lightest supersymmetric particle (LSP). 

-T— , , We consider the unique hadron collider signatures of a weak scale right-handed sneutrino LSP, 

H assuming i?-parity conservation. For concreteness, we concentrate on stop pair-production and 

decay at the Tevatron and the Large Hadron Collider, and briefly comment on the production and 
decay of other supersymmetric particles. 
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I. INTRODUCTION 



A renormalizable extension of the standard model (SM) that incorporates neutrino 
masses can be obtained by introducing (at least two) SM gauge singlet fermions - right- 
handed neutrinos Nji. Neutrino masses arise due to the fact that right-handed neutrinos 
can couple to the left-handed ones (z/^) through Yukawa couplings. After electroweak 
symmetry breaking, the SM neutrinos are endowed with either Majorana or Dirac masses. 
A supersymmetric version of this scenario implies the existence of new, complex, SM gauge 
singlet scalar fields - the right-handed* sneutrinos Nji. 

The right-handed neutrino superfield iV is a SM gauge singlet and is allowed to have 
a supersymmetry (SUSY) preserving Majorana mass M^. In the usual high-scale seesaw 
mechanism Mn - the scale of lepton number breaking - is, in general, not related to the 
SUSY breaking scale. It is often considered to be around 10^'* GeV so that, if the neutrino 
Yukawa couplings are order one, active neutrino masses are around 0.1 eV, the scale inferred 
from neutrino oscillation experiments 2[. This scenario, while elegant, is very hard to verify 
experimentally^ and is motivated by the fact that neutrino Yukawa couplings are "naturally" 
expected to be of order one. If the seesaw energy scale is indeed very high, right-handed 
sneutrino masses are expected to be of order Mn and, hence, these are decoupled from low- 
energy phenomena. SUSY effects in a theory with a high seesaw scale have been considered 
in ^. It is also possible that neutrino masses arize from couplings to an SU(2) triplet 
Higgs boson [5]. The supersymmetric version of such a theory doesn't have to include a 
right-handed neutrino superfield and thus its signatures is outside the purview of this work. 

It has recently been emphasized that any value of is technically natural Gj (including 
Mtv = 0, in which case the neutrinos are Dirac fermions), and that it is important to 
explore the phenomeno logical consequences of all values of Mtv. One intriguing possibility 
is to consider, for example, that the unknown physics of SUSY breaking and lepton number 



* "Sfermions" are, of course, scalar particles, and have no sense of handedness. Throughout, however, we 
refer to left- and right-handed sfermions (scalar neutrinos, scalar tops, etc), as is commonly done in the 
literature, in order to indicate the super-partner of the various left- and right-handed chiral fcrmion fields 
(the neutrino, the top quark, etc). 

One possibility is to look for right-handed neutrino traces in the RGE evolution of the Soft SUSY breaking 
parameters, which could be revealed using precision measurements at a next-generation linear collider [3|. 
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breaking are intimately connected, resulting in a common mass scale. If this is the case (or if 
the SUSY breaking scale were larger than Mn), right-handed sneutrino masses are expected 
to be of the same order as the masses of all other super-partners. Finally, the radiative 
stability of the weak scale indicates that the SUSY breaking scale ought to be around the 
weak scale, so that super-partners are expected to make their presence felt at TeV scale 
collider experiments. 

Depending on the physics of SUSY breaking, right-handed sneutrino masses may end 
up below all other "active" super-partner masses, so that the lightest sneutrino is the 
lightest supersymmetric partner (LSP). In general, left- and right-handed sneutrinos mix. 
However, given that the seesaw energy scale is assumed to be very low, the neutrino Yukawa 
coupling is, phenomenologically, required to be very small. Hence, if one assumes that the 
supersymmetric breaking "A-terms" are proportional to the respective Yukawa couplings, 
left-right sneutrino mixing is expected to be very small (see, for example, for details, 
some aspects of which are given in App. [A]for convenience) and the LSP will turn out to be 
composed of a mostly right-handed sneutrino, Nn. 

In a recent work Q], we considered such a Nr LSP when the Majorana and SUSY breaking 
mass scales are around the electroweak scale. In this case, the neutrino Yukawa coupling 
Yn ~ 10"^ in order for the light neutrino masses to be to be of the order of 0.1 eV. We 
analyzed the sneutrino sector, pointed out in which cases the left-right sneutrino mixing angle 
can be tiny resulting in an almost pure LSP, and argued that the Nn is an interesting 
non-thermal dark matter candidate^ (see also [3]). In this paper we explore the hadron 
collider implications of a predominantly Nfj^ LSP. Collider signatures of a mixed ul — Nn 
LSP have been explored previously in other studies [lo| . 

Since a weak scale Nr LSP interacts only through the tiny Yukawa coupling Yn, many 
unique collider signatures are expected to arise. We will show that the most striking among 
them is the possibility of the next-to-lightest supersymmetric particle (NLSP) being long- 
lived enough to leave a displaced vertex in the detector. A similar situation exists in theories 



of gauge mediation where the gravitino is the LSP UJ, with usually the scalar tau or 



the lightest neutralino as NLSPs. A displaced vertex signature in the context of "hidden 



^ In case of additional non-standard interactions, it can also be a thermal dark matter candidate 



2 



valley" models is discussed in [l3|. However, the scenario discussed here is distinguished by 
(a) the potentially long-lived states can be strongly interacting, (b) the LSP carries lepton 
number, which implies associated leptons in the final state, and (c) the LSP interacts only 
through the Yukawa coupling Y^^ which generically implies non-universal rates for e, /i, and 
T type leptons. 

In order to illustrate the unique collider signatures of a Nr LSP, we consider in detail 
the case where the right-handed stop (t/j) is the NLSP^, and analyze its pair-production 
and decay. While serving to illustrate many of the unique features of a Nr LSP, a light 
stop is natural in many scenarios {14}] and is favored in successful theories of electroweak 
baryogenesis [l5|. Stops, being strongly interacting particles, are also expected to be 
produced at significant rates at the Tevatron and the LHC Stop production and decay 
have also been analyzed in other contexts 111, Q, [l7 |. 

The outline of the paper is as follows: In Sec. [Tll we compute the decay rate of the into 
the 3-body final state M'^Nr- We incorporate this decay matrix element into the Monte 
Carlo program Pythia (version 6.327) ISj and study the Tevatron and LHC signatures 
in Sec. IIIII In Sec. HVl we briefiy discuss the dominant signatures of other SUSY NLSP 
candidates, such as gluinos, sbottoms, gauginos, and sleptons, and comment on "co-LSP" 
right-handed sneutrinos. We offer our conclusions in Sec. |Vl The details of the model with 
which we work were spelled out in detail in [3], and are summarized in App. |X1 The exact 
expression for the stop decay matrix element is given in App. [Bl 



II. tR DECAY TO Nr 

In the scenario of interest, assuming i?-parity conservation, all supersymmetric particles 
eventually cascade-decay to the stable Nr_ LSP. Since the Nr couples only via the tiny Y/v, 
the NLSP can be potentially long-lived. To illustrate this aspect, in this section we consider 
in detail the decay of the right-handed stop {Ir). This will be particularly relevant when 
the Ir is not too heavy, in which case its production cross-section can be large enough to 

^ For simplicity we will take the light i to be predominantly iji. While stop mixing is, in general, not 
negligible, we find that its inclusion does not change, qualitatively, any of our results. 
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FIG. 1: tR biNji decay mode. The arrows indicate fermion number flow. 



be observable at the Tevatron or the LHC A brief discussion of the decays of other NLSP 
candidates is given in Sec. IIV[ We omit, for simphcity, stop mixing. We find that the 
inclusion of left-right scalar top mixing does not change significantly any of our results. 

Fig. [1] shows the dominant contribution to the decay of a predominantly right-handed 
stop, which is assumed to be the NLSP, into a pure Nji LSP. The related mode t —>■ tuNji 
can also be relevant if kinematically allowed, although phase-space suppressed with respect 
to t — i> MNr due to the top quark in the final state. 

The complete expression for the lightest stop decay matrix element, including the mixings 
of all relevant SUSY particles, is presented in App. [Bl For a purely right-handed stop and 
sneutrino, in the limit of pure higgsino H exchange^, the formula for the stop decay width 
has a simple form. In the t/j rest frame, defining 6be as the angle between the bottom quark 
and the charged lepton, the matrix element squared is given by 



Here, \Yt\ {\Yn\) is the top (neutrino) Yukawa coupling, (Ei) is the 6-quark (charged 
lepton) energy, while kb is the 6-quark four-momentum. Mj^ and pi^ are, respectively, 
the ifi mass and four-momentum, while is the Higgsino mass. Due to the (1 + cos^;,^) 
factor, the matrix element peaks when the 6-quark and the charged lepton are aligned. Later 
(Sec. IIII|) . we will point out that this cos 6bi behavior can be used to distinguish between 
tfi — > MNfi and t bW —>■ Mu. It turns out that, for top quark decays, there is no peak in 
the event distribution when the 6-quark and the charged lepton are aligned (c/Fig. [7]). 

^ We find that, for the intentions of this paper, it is safe to ignore higgsino-gaugino mixing. 




A\Yt\^\Y,,\^MlEbEe {1 + cos9m) 



(1) 
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FIG. 2: SOLID - Contours of constant cr for the right-handed stop, in the Mj^ — plane. We 
assume Mp^ = l.lMj^. The Mj^ = ^fj^ dotted hue bounds the region of parameter space where 
is the LSP. DASHED - 3 u sensitivity for stop pair production at the Tevatron (left) and the 
LHC (right), for different integrated luminosities. Note that the displaced vertices are not taken 
into account when defining the sensitivity. See text for the details. 



The decay rate is given by 

.Ml 

' .(16vr2) 

where f^ps is a dimensionless 3-body phase space function. Note that > Mi^. 



r = 4|yj2|F^|2— ^ 



47r ^ 

72 hps 



(2) 



Numerically, for = 100 GeV, 



0.05 



„ . 10 mm . 225 GeVr /^^^V 

We remind the reader that f^ps is a function of M^^, M^^, and M^. 

For Yi\i ~ 10~^, cr can vary between a few millimeters to several meters, depending on 
the stop, sneutrino, and Higgsino masses, cr values in the M^^ — Mfj^ plane are depicted in 
Fig. [21 assuming that Mfj = I.IM^^. For larger values of the Higgsino mass, the constant cr 
contours move toward the Mj^-axis. Displaced vertices, z.e., cr ^ 1 mm are to be expected 
even for large scalar top masses, while for light enough stops and heavy enough LSPs and 
Higgsinos the stop may even be collider-stable. We elaborate on these possibilities in the 
next section. 

If lepton number is conserved in nature, = = ce = (see App. neutrinos are 
Dirac fermions and, in order to obtain the right order of magnitude for the active neutrino 




/ 
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FIG. 3: Dominant parton level t pair-production processes at a hadron collider. The diagrams 
with gluons in the initial state dominate at the LHC. 



masses, ~ 10 values are required (c/ Eq. flA2p ). In this case, cr is measured in 
kilometers and, as far as collider phenomenology is concerned, the t/j is absolutely stable. 

Heavy, strongly interacting, collider-stable SUSY particles are expected to form 'R- 
hadrons,' which behave like heavy nucleon-like objects. Similar experimental signatures 



have been studied elsewhere 



19| . We do not pursue other experimental signatures of heavy. 



long-lived, had^nc state. We would, however, Uke to co„.„.ent on the fact that vety long- 
lived tfi could form narrow "onium-like" t/jt^ bound states** [211], whose decay may lead to 
a peak above the continuum background. An investigation of this phenomenon will be left 
for future work. 



III. tR AT THE TEVATRON AND LHC 

We consider the production and decay of a tji NLSP at the Tevatron and the LHC followed 
by its subsequent decay into a purely right-handed sneutrino N^i. LSP, i.e., — > M^Nr. 
Unless noted otherwise, we assume that the Ir decays in this way 100% of the time. To 
illustrate the behavior of various observables, we will consider the case Mi^ = 225 GeV, 
M^^ = 100 GeV, = 250 GeV, 1^ = 4 x 10"^. In this case, cr ^ 10 mm (see Eq. ([3])). 

Fig. [3] depicts the dominant production mechanism of a Ir pair*, at leading order. The 
** The same would also be true of long-lived, gluino bound states [2o| . 

* Depending on the details of the SUSY spectrum, there may be other tji production channels, including 



6 



TABLE I: "Level 1" cuts imposed in our analysis. We define R^^ = {(pb — + {r]b — Ve)'^- 



1. 


Rapidity cuts 


Irjil < 2.5, < 2.5 


2. 


Pt cuts 


PTi > 20 GeV, pTb > 10 GeV 


3. 


Isolation cut 


Rbe > 0.4 



stop Disp Vtx (ctau = 10 mm) stop Disp Vtx (ctau = 10 mm) 
0.14 I 1 1 , 1 1 0.14 I 1 , 1 




dTvx (mm) dTvx (mm) 

FIG. 4: The distributions of the transverse displacement of the stop (in mm), at the Tevatron 
(left) and the LHC (right). 



dominant contribution to t/j decay is depicted in Fig. [H The escapes the detector 
as missing energy so that, experimentally, one should observe pp{p) — > ^ijt/j h[^M~ + 
$j, . The dominant physics background for this process is expected to be top quark pair 
production: pp{p) ^ ti ^ hW^hW~ bi~^M~ + , where the missing energy is due to 
neutrinos in the final state. We use the Monte Carlo program Pythia (version 6.327) in order 
to analyze, at the parton level, the signal and background at the Tevatron and the LHC. 

The observables available in order to discriminate the signal from background are the 
transverse components of k^, k^+, k^ and k^-. We impose the "level 1" cuts shown in Tabled 
necessitated by the detector geometric acceptance and thresholds. For the isolation cut, we 
make the standard definition Rl^ = {(pb — (p^Y + {rjb — rji)'^. These values are typical of other 
Tevatron analyses, and are meant to approximate the capabilities of the LHC detectors. 

gluino pair production, followed by g ^ ttji (if kinematically accessible) . Here we concentrate on "direct" 
QCD production. 
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Fig. m depicts the resulting transverse displacement (in the x—y plane), after including the 
boost of the stop at the Tevatron (left) and the LHC (right). Even though the stop is more 
boosted at the LHC (compared to the Tevatron) due to the higher center-of-mass energy, 
the stop pair is more forward peaked at the LHC, causing their transverse displacement to 
be smaller at the LHC (compared to the Tevatron). For other and Mi^ values, the 

transverse displacement scales with the scalar top's cr. These have been discussed in the 
previous section, and are depicted in Fig. [2l 

If a stop is long-lived enough, it will behave like a stable or quasi-stable hadronic object (as 
commented on earlier). If it decays before exiting the tracking subsystem, a displaced vertex 
may be reconstructed through the stop decay products' 3-momenta meeting away from the 
primary interaction point. In the example considered, a displaced vertex of 10 mm can be 
easily discerned at the Tevatron and the LHC. On each side, the 6-quark itself leads to an 
additional displaced vertex, and its 3-momentum vector can be reconstructed from its decay 
products. In combination with the 3-momentum of the lepton, the stop displaced vertex can 
be determined. In order to reveal the displaced vertex, one must require either the 6-quark 
or the charged lepton 3-momentum vector to miss the primary vertex. Since a pair of stops 
is produced, we would expect to discern two displaced vertices in the event (not counting 
the displaced vertices due to the b-quarks). If a 6-quark (£+) cannot be distinguished from 
a 6-antiquark (^~), the associated combinatoric problem of assigning the decay products to 
the parent stop will have to be considered. Such an event with two displaced vertices, from 
each of which originates a high pt C- and 6-quark, is quite uncommon in SUSY models^ and 
might prove to be one of the main distinguishing characteristics of such a scenario. A cut 
on the displaced vertex provides for a very effective way to separate stop events from the 
top background. If one can efficiently explore such cuts, we anticipate that NLSP scalar top 
searches may turn out to be physics-background free. 

Another characteristic feature is the non-universal rates for decays into e, /i and r leptons. 
This is to be expected given that the stop decays are proportional to the Ya^'s which are in 
general different for the three leptons. Here, we do not take advantage of this feature. 



'I' Other SUSY models that lead to a displaced vertex are the NMSSM with a singlino LSP [25(, the case 



of bilinear i?- parity breaking where neutrino masses are generated via i?-parity breaking effects 26| , and 



also those in 27 1 
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FIG. 5: The distribution of px of the 6-quark at the Tevatron (left) and LHC (right) resulting 
from the decay of a 225 GeV stop and a top quark. 

If the stop displaced vertex cannot be efficiently resolved, one will have to resort to 
more conventional analysis methods. In the remainder of this section, we explore various 
kinematical distributions for both the signal (right-handed scalar top pair production) and 
the physics background (top pair production), obtained after imposing the level 1 cuts listed 
in Table [H All the distributions are in the lab frame, normalized to unit area. Note that the 
analysis performed here also applies to other SUSY scenarios in which the scalar top decays 
predominantly to bottom plus charged lepton plus missing energy, regardless of whether the 
stop decays promptly or leaves behind a displaced vertex. 

Fig. [5] depicts the parton level distributions of the pt of the 6-quark at the Tevatron 
(left) and the LHC (right) resulting from a 225 GeV stop and from a top, obtained using 
Pythia. Fig. E] depicts the distributions of the pt of the charged lepton at the Tevatron 
(left) and the LHC (right), resulting from a 225 GeV stop, and from a top. The pt of 
the 6-quark from the 225 GeV stop peaks at a lower value compared to the top quark 
background, and therefore accepting them at high efficiency for p^ ^ 40 GeV will be very 
helpful in maximizing the signal acceptance. The signal and background shapes are quite 
similar and no simple set of pt cuts can be made in order to significantly separate signal from 
background. As the mass difference between the and the Nr decreases, the 6-quark and 
the charged lepton px distributions peak at lower values due to less available phase-space. 
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cos(th)_bE cos(th)_bE 




I ' ' ' ' ' ' ' ' ' ^ 1 1 1 1 1 1 1 1 1 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cos(th) cos(th) 



FIG. 7: The distribution of cosOm resulting from the decay of a 225 GeV stop, and a top, at the 
Tevatron (left) and LHC (right). 



making the measurements more challenging. 

Fig. [7| depicts the distribution of cos^^^, the angle between the 3-momenta k;, and k^, 
for both the signal and background. It is important to appreciate that, by default, Pythia 
generates stop decays into the 3-body final state according only to phase-space, ignoring 
the angular dependence of the decay matrix element. We have reweighted Pythia events to 
include the correct angular dependence in the decay matrix element. Consistent with the 
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FIG. 8: The total jy^ distribution resulting from a 225 GeV stop, and a top, at the Tevatron (left) 
and the LHC (right). 



expectation from Eq. ([T]), we see for the signal that the distribution peaks for the 6-quark 
and charged lepton 3-momenta aligned, unlike the background. -f It is unfortunate that the 
isolation level 1 cut (see Table [T]) on the leptons removes more signal events than background 
events. Relaxing this constraint as much as practical would help in this regard. 

Fig. [8] shows the total distribution resulting from the production and decay of a 
225 GeV stop and a top quark for the Tevatron (left) and the LHC (right). Here />p is 
defined as \//>r ^ + /"r ^, with prp ^ and j^^^ the x and y components of the total missing 
momentum vector. Pythia ignores the spin correlation between the two opposite side 
particles and therefore the />p distributions shown are not entirely accurate in case of the 



t-quarks. However, spin correlation modifies the distribution to only a small degree 22| , 
although it can lead to significant effects in suitably chosen observables (23|. 



The angular correlation of the stop pair is different from that of the top quark pair, since 
the former is a scalar and the latter a fermion. We expect the quantities kf, ■ and k£+ ■ k^- 
to inherit some of this difference, making these potentially good discriminants for signal and 
background. As just mentioned, Pythia is not suitable to investigate this aspect since it 
does not retain spin correlations, and we postpone this investigation to future work. 



For another discussion on how to extract stop NLSP's (in a different context) see 11 1 
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TABLE II: Signal (stop) and background (top) pair production cross-section at the Tevatron, with 
eb = 0.5, and = 0.9. as and as denote the signal and background cross-sections, a the fraction 
that passes level 1 cuts (see Table |I|), S and B the number of signal and background events for 
1 fb~^ 







0-5 (pb) 


o-B(pb) 


a 


S 


B 


S/B 


s/Vb 


S/^/S + B 


100 


50 


11.83 


6.77 


0.26 


162 


9 


18.93 


55.36 


12.4 


75 


11.83 


6.77 


0.04 


4 


9 


0.45 


1.31 


1.09 


150 


100 


1.24 


6.77 


0.29 


21 


9 


2.46 


7.21 


3.87 


125 


1.24 


6.77 


0.05 


1 


9 


0.07 


0.21 


0.21 


175 


100 


0.48 


6.77 


0.47 


22 


9 


2.53 


7.39 


3.93 


150 


0.48 


6.77 


0.05 


0.2 


9 


0.03 


0.08 


0.08 


250 


100 


0.04 


6.77 


0.71 


4 


9 


0.48 


1.4 


1.15 


200 


0.04 


6.77 


0.31 


1 


9 


0.09 


0.27 


0.26 



TABLE III: Signal (stop) and background (top) pair production cross-section at the LHC, with 
eb = 0.5, and ee = 0.9. as and as denote the signal and background cross-sections, a the fraction 
that passes level 1 cuts (see Table Hj), S and B the number of signal and background events for 
10 fb'\ 







crs(pb) 


o-B(pb) 


a 


S 


B 


S/B 


S/^/B 


S/^S + B 


100 


75 


1332.5 


873 


0.03 


1938 


8662 


0.22 


20.82 


18.82 


83 


1332.5 


873 


0.01 


73 


8662 


0.01 


0.78 


0.78 


150 


100 


228.79 


873 


0.23 


25325 


8662 


2.92 


272.1 


137.37 


128 


228.79 


873 


0.02 


144 


8662 


0.02 


1.54 


1.53 


250 


200 


21.32 


873 


0.26 


2886 


8662 


0.33 


31.01 


26.86 


225 


21.32 


873 


0.03 


40 


8662 


0.01 


0.43 


0.43 


500 


400 


0.56 


873 


0.59 


398 


8662 


0.05 


4.28 


4.19 


425 


0.56 


873 


0.48 


263 


8662 


0.03 


2.83 


2.78 


650 


250 


0.14 


873 


0.83 


195 


8662 


0.02 


2.10 


2.08 


500 


0.14 


873 


0.71 


145 


8662 


0.02 


1.56 


1.54 
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The sensitivity of the Tevatron and the LHC to the stop NLSP is depicted in Table |TT] 
and Table llllt respectively, for various and Nr masses. The stop pair-production cross- 



section from Pythia multiplied by the appropriate K-factor 2j] is displayed in the table 



for the Tevatron and the LHC, along with the top pair-production cross-section [28j. The 
fraction of signal events that passes level 1 cuts (specified in Table [T]), denoted a, obtained 
using Pythia is shown, and the fraction of background events (not shown in the table) is 0.79 
(Tevatron), and 0.7 (LHC). We compute the number of signal events (S*), background events 
(5), for 1 fb~^ (Tevatron) and 10 fb~^ (LHC), and compute the figures-of-merit S/ B, S/\fB 
and S/ \/S + B. We have taken the b-tagging efficiency = 0.5 and the lepton identification 
efficiency = 0.9. The number of events is given by crCa'^elej, where C is the integrated 
luminosity. It is important to emphasize that we do not consider the possibility that one 
can efficiently identify that the scalar top decays far from the production point. If this is 
the case, we anticipate that physics backgrounds will be significantly reduced. At the same 
time, we remind the reader that we have fixed Mfj/M^^ = 1.1. For larger values of this 
ratio, scalar tops are expected to be (much) longer-lived, in which case other approaches to 
data analysis are required. 

Fig. [2] shows the projected 3 a contours in the Mj^ — plane, with 1 and 10 fb~^ at 
the Tevatron^, and 10 and 500 fb~^ at the LHC. We find from our analysis with just level 1 
cuts (see Table [T]) that the Tevatron can probe stop masses up to about 300 GeV (with 
10 fb~^), while the LHC reach extends slightly above 650 GeV (with 500 fb^^), with the 
reach depending on the Nji mass. We expect that the sensitivity can be improved with more 
sophisticated cuts and a multivariate analysis. A smaller M^^ — mass difference leads to 
softer 6-quarks and charged leptons, resulting in a smaller a and leading to a lower statistical 
significance. Here we assume that the stop decays to the leptonic channel considered with 
100 % branching ratio. In a specific model the actual significance can be obtained by 
including the branching ratio. Our estimates for the top and stop production rates are 
in good qualitative agreement with experimental measurements of top production [3.c| and 
bounds imposed by stop searches at the Tevatron [si^. 



The i — > Mi>L mode at the Tevatron has been studied in the j£^£ l^^ channel in [29] with quahtatively 
similar results to ours. 
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IV. PRODUCTION AND DECAY OF OTHER SUSY PARTICLES 



In this section we briefly discuss some aspects of unique signatures associated witfi 
tlie production and decay of otlier SUSY particles. The actual signatures depend on the 
detailed SUSY spectrum. In particular, we digress on features associated with other NLSP 
candidates, including sbottoms, gluinos, gauginos and sleptons. We also comment on the 
fate of the "other" mostly right-handed sneutrinos, which are expected to be, as far as 
collider experiments are concerned, also stable. 

A. sbottom 

If the is the NLSP, its pair-production cross-section is naively similar to the stop 
NLSP case discussed in Sec. IIIII The decay channel analogous to the one discussed for the 
stop NLSP is 6 ^ tiNn. This decay channel, however, may not be kinematically accessible. 
Regardless, the dominant channel for a large chunk of the parameter space is expected to 
be be bu buN^, resulting in the signature pp{p) bb + $rp . Due to the additional 
suppression of the decay rate by Y^, we expect NLSP sbottoms to be more long-lived than 
the stop NLSP by a factor of /Y^. Here, we expect that it will be harder (compared to 
the NLSP case) to identify efficiently a potentially displaced vertex, due to the absence 
of charged leptons in the final state. One may be able to achieve this by asking whether the 
reconstructed 3-momenta of the opposite side 6-quarks point back to the primary vertex. One 
source of background is QCD bb production - huge - but demanding substantial should 
help extract the signal. We expect that substantial Ipj^ will also be crucial for triggering on 
the event. 

The decay mode 6/? — > clN^ is suppressed by \Vcb\^ compared to the mode discussed 
above, and should have a branching ratio around 10~^. This mode leads to the signature 
PP{p) ~^ cci'^i~ + ]^rp . The additional pair of leptons will help in discriminating these events 
from background but, given the small branching ratio, we expect this decay mode to be out 
of the reach of the Tevatron (but not the LHC). As with stop NLSPs, given that the decay 
is proportional to the Yat's, we generically expect different rates for different charged lepton 
final states as a characteristic of this scenario. 
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B. Gluino 



If the gluino is the NLSP, it decays primarily to a four-body final state g 



qq 



Nr, via 



off-shell q' and x^. This leads to the signature pp{p) —>■ 4:j + 21"^ + ^j. . A related mode is to 
a neutrino in the final state, i.e., g — > qquNji via off-shell q and x^, leading to the signature 
PP{p) ~^ 4j + Jprp . Unlike the previous mode, in this mode the leptons are unobservable. 

Owing to the four-body final state and the Yat suppression, the gluino is expected to be 
quite long lived. The decay rate is 

Air 



2 2-1^2 g 
9s9 ^Nj^2M^ 



(levr^ 



Zsf^PS 



(4) 



where f^ps is a dimensionless four-body phase-space function. Compared to the NLSP 
decay width (Eq. [2]), 



CTa 



3PS 



(5) 



4P5. 



We naively estimate Tg ^ so that, for all practical collider purposes, the gluino 

is stable (cr > 100 m). Such a long lived gluino forms an R-hadron and some of its 
experimental signatures have been discussed in 1^. 

As an aside, we comment that if the mass spectrum is such that the gluino is no longer 
so long lived, the fact that the gluino is a Majorana particle can be used in order to identify 
gluino production and decay. Such processes can result in decays into a pair of same-sign 
leptons, and, if kinematically accessible, the case of decays into same-sign tops can be easily 
distinguished from background [17 1. 



C. Gaugino 

Charginos can be pair-produced via off-shell 7 and Z-boson exchange, and the NLSP 
chargino decays into the LSP via x^ —>■ I'^Nr. This leads to the signature pp{p) t^i'^ + 

. The production cross-section is suppressed relative to that of strongly interacting SUSY 
particle pair production by {g/ gsY ~ 10~^, and therefore probably too small to be probed 
at the Tevatron. The LHC, on the other hand, will have the ability to produce weakly 
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interacting states in significant numbers. The chargino lifetime can be easily estimated as 

(6) 

where f2ps is a dimensionless two-body phase-space function. Comparing this with the stop 
NLSP decay width (Eq. [2]), we obtain 

For Mfj = I.IM^^, f2ps ~ 0(1), and our numerical estimate f^ps ~ 0.05, we get cr^ ~ 
cr^-^ X 10~^. Chargino decay is, as far as hadron collider experiments are concerned, prompt. 
The main background is expected to be H^-boson pair production. 

Neutralinos can be pair-produced via an off-shell squark, and an NLSP neutralino decays, 
most of the time, invisibly: vNji. A gluon jet or photon can be radiated from the initial 

state leading to the signature pp{p) i$x or 7^^. . Due to the electroweak production 
cross-section, it is unlikely that such processes are accessible at the Tevatron. At the LHC, 
on the other hand, one may run into a large number of these events. Whether or not 
these can be extracted from the various backgrounds requires a dedicated study, beyond the 
ambitions of this brief section. 



4:71 

167r^ 



2PS 



D. Slepton 

Sleptons are pair produced by the exchange of virtual 7 and Z-boson exchange, followed 
by the decay £ ^i^Nji. This leads to the signature pp{p) — > + $x . Since this is a 
three-body decay, like in the case of the stop NLSP, we expect displaced vertices or very 
long-lived sleptons. The only observable decay product of the NLSP slepton is the charged 
lepton, so the displaced vertex is characterized by a lepton track that does not point back to 
the primary vertex. The dominant background for this channel is VT-boson pair production. 
Note that this case is similar to some manifestations of gauge-mediated SUSY breaking with 
the gravitino as the LSP and the slepton (usually the scalar tau), as the NLSP |l2| . 
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E. Co-LSP right-handed sneutrino 



Given that there are at least two generations of right-handed sneutrinos (N^^), we explore 
the possibility of observable consequences of these "co-LSP's." Such a possibility was already 
raised during studies of the lightest left-handed sneutrino as the LSP [32]. It is possible, 

~ (2) 

in all the decays considered so far, that a heavier co-LSP, say N}^ , is produced. This 
state later decays to the "real" LSP, say The relevant (observable) decay channels are 
N^^ — > N^^i~^i~ and the one-loop decay N^^ ^r^I- An estimate of these decay widths 
is 



rv(2) ~ J AT , ^9 9\9 i3P5' 




(8) 



^ /3P5GeV 



where fsps is a dimensionless phase space factor. Note that this decay width is proportional 
to Yn - tiny! - to the fourth power. This leads to a lifetime well above 1 s, so that N}^ is 
long-lived enough to exit the detector unseen. 

A concern for such long lived particles is whether they disrupt the successful predictions 
of big-bang nucleosynthesis. Note that, if the mass difference between the co-LSP states is 
relatively small (as one would naively expect), fsps ^ 1, and the lifetime may be orders 
of magnitude longer than the naive estimate above. A more more detailed study, beyond 
the ambitions of this section, is required in order to determine the impact of these decaying 
co-LSPs in the early universe. 



V. CONCLUSIONS 

Right-handed sneutrinos [Nji) are present in a supersymmetric theory that includes right- 
handed neutrinos. If the right-handed neutrino Majorana mass and the SUSY breaking 
scale are identified with the electroweak scale, our understanding of light neutrino masses 
requires the neutrino Yukawa coupling to be ^ 10^^. Depending on the details of the 
SUSY breaking mechanism, it is possible that a predominantly right-handed sneutrino is 
the LSP. We detailed such a theory in Q], and explored the cosmological implications of 
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a right-handed sneutrino LSP. In this paper we study some hadron coUider (Tevatron and 
LHC) signatures of such a Nr LSP. 

We show that if such a Nji is the hghtest supersymmetric particle (LSP), the colhder 
signatures are very interesting, owing to the fact that the Nji interacts only through the 
tiny Y/v- If -R-parity is conserved, all heavier SUSY particles eventually have to cascade 
decay to the Nn through the interaction parameterized by Y^. Among other potentially 
observable effects (leptons, missing energy), we find that the next-to lightest supersymmetric 
particle (NLSP) is potentially (very) long-lived. If a mostly right-handed sneutrino is the 
LSP, one generically expects displaced vertices due to the relatively long NLSP lifetime, or 
heavy, collider-stable hadronic or weak states. Since the NLSP decay to the sneutrino LSP 
proceeds through a Yukawa coupling, we expect non-universal rates in the e, fi and r lepton 
channels. Even if the right-handed sneutrinos constituted all of the dark-matter observed, 
given its tiny interaction cross-section, the event rate in direct-detection experiments will 
be well below the sensitivity of planned experiments. These aspects help in distinguishing 
this scenario from other SUSY models in interpreting any deviations from SM predictions. 
Of course, in order to decide whether the scenario discussed here is realized in nature, it 
would be imperative to establish that low-energy SUSY is indeed realized in nature. In 
order to achieve this, it is important to observe a handful of superpartners, and establish 
that their spins and interactions agree with the predictions of a supersymmetric version of 
the standard model. A detailed discussion of this is beyond the ambitions of this paper, but 
we refer readers to |33| for a list of very recent discussions of this issue. 

In order to illustrate these aspects, we consider, in Sec. [TTl the scalar top 3-body decay 
tfi — i> M~^Nji via an off- shell H. We calculate the decay matrix element showing explicitly 
the dependence on 6^, the angle between the 3-momenta of the final-state 6-quark and 
charged lepton. We show that for ~ 10~^, cr for the scalar top varies from millimeters 
to meters. 

We performed a simulation of stop pair production and decay at the Tevatron and the 
LHC, using the Monte Carlo program Pythia. Pythia by default treats the 3-body final state 
according only to phase-space, so we modified the program to correctly incorporate the decay 
matrix element. The signature of this signal is pp{p) tRt*R M^hl" + . The dominant 
physics background is top quark pair production pp{p) ^ tt ^ bW'^bW^ M^M' + $rp , 
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where the missing energy is due to final state neutrinos. 

The transverse displaced vertex after folding in the boost of the stop is shown in Fig. HI 
and is easily discernible at the Tevatron and the LHC This can be used to very effectively 
suppress the background. If the parameters are such that a displaced vertex cannot be 
resolved, we will have to rely on the distributions and event rates in order to show an excess 
above background. We see that the px distributions of the 6-quark and the charged lepton 
from a 225 GeV stop and the top background are quite similar, and no clear set of pt cuts 
can be applied to separate them. The cos 6m distribution is, however, sufficiently different 
for signal and background, as can be seen in Fig. [71 

We summarize the Tevatron and LHC reach in Tables [III and IIIII From this physics level 
study, we estimate that the Tevatron can probe stop masses up to about 300 GeV with 
10 fb~^, while the LHC is sensitive to about 650 GeV with 500 fb^^ of integrated luminosity. 
This reach was obtained after applying only the level 1 cuts shown in Table HI A more 
sophisticated analysis is expected to extend the estimates obtained here. Needless to say, 
a full detector-level simulation is necessary in order to realistically assess the capabilities 
of the Tevatron and the LHC. We would like to point out that the results presented in 
Tables [TTl and IIIII should also apply to other SUSY scenarios in which the scalar top decays 
predominantly to a 6-jet, a charged lepton, and missing transverse energy, regardless of 
whether the stop decay occurs promptly or leaves behind a displaced vertex. 

Of course, the nature of the NLSP depends on the details of supersymmetry breaking. 
We offered some remarks in Sec. [IVl on what channels are promising for several potential 
NLSP candidates. We leave a comprehensive analysis of many possible cascade decay chains 
for future work. 
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APPENDIX A: THE MODEL 



To the field content of the MSSM, we add (for each generation) a right-handed neutrino 
superfield = {Nr, N, Fj^). Written as left-chiral fields, the superfields are: Q, U^, D^, L, 
E^, N'^. As usual, the MSSM Higgs doublet superfields are IL. and Hd- Here we repeat the 
main aspects of the model, details of which can be found in [7|. 

The superpotential is 

W = Lf^YuQ ■ Hu - D'YnQ ■ Ha + N'YnL ■ % - E'YeL ■ % + N'^N' + iiH^ ■ , (Al) 

where A ■ B denotes the antisymmetric product of the fields A and 5, Y are the Yukawa 
couplings that are 3x3 matrices in generation space, and, M^v breaks lepton number. 

After electroweak symmetry breaking (when the Higgs scalars get vacuum expectation 
values Vu and Vd), the lowest neutrino mass eigenvalue is given by the standard seesaw 
relation 

= ^ , (A2) 

where we assume VuYn <^ M^r. Neutrino oscillation experiments indicate that niy ~ 0.1 eV, 
and if ~ Eq. ([M]) implies that Y^ ^10"^. 



The soft SUSY breaking Lagrangian is given by 

C-susYBr = - qi^lqL - u'^mluR - d^j^mldR + (-uJjA^gi ■ K + d^^AdqL ■ hd + h.c.) 

- i{m]iL - N^mlNR - e^RmlcR + (-iV^A^4 ■ K + 4^e4 ■ hd + h.c.) 
+ Vl-h^Y'l{l.K) + Nt-^NR + h.c. 

+ (6/x/i„ ■ hd + h.c.) , (A3) 



where q and b^M^ are SUSY-breaking, lepton-number breaking parameters. As usual, 
are SUSY breaking scalar masses- squared, A are SUSY breaking A-terms, and b is the SUSY 
breaking Higgs boson B-term. 
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After electroweak symmetry breaking, the sneutrino mass matrix (generation structure 
suppressed) is given by 



/ 



where mli^ 



mil 



-vice 
\vuMI,Yn 



m 



2t 
RL 



^RR 
-hNMN 



vuMj;y* 



VuY^Mm \ 



m 



m 



2* 
LL 

2 * 
RL 



m 



m 



2T 
RL 

2* 
RR 



J 



m 



2 

RR 



(A4) 



(MjvM^ + m% + vIYnYj^), m|i = (-^*t;dFjv + 
VuA^), and = (m|/2) cos 2/5 is the D-term contribution. 

The vl-Nr mixing angle is given by (see [3] for details^) 

_ 2 (m|^ ± VuM^^Yn) 

{^LL T vice) - [ml^ T hNM^) 

If Atv oc Yat, as is the case in several popular SUSY breaking scenarios, it is easy to see that 
sin^p ~ Yn. For Yn ~ 10~^ - the case of interest here - the mixing angle is tiny, and, as 
long as < the LSP is an almost pure Nr. 



APPENDIX B: THE DECAY ii hl+Di 

Here, we present formulas for the decay width of the 3-body decay of the stop, ti — > hl^vi. 
Although in the main body of the paper we consider the case when ^ 1 with z/i Nr 
as the LSP, the formulas in this section are valid in general. 

The matrix element is given by 

Tf, = Y: vM (kliPL + 1\iPr) -r-^ {KiPr + C^l) v.{p,) , (Bl) 

where 

l\i = -9 cos e^Vii - Yt sin eiVi2 , /c-i = Yfe cos eiUi2 , 

(B2j 

= -gcosOyVii - Y^sm.6yVi2 , fc^^ = Y^ cos 6'£;[/j2 , 

in the case of Dirac neutrinos. In the case of Majorana neutrinos, l^^ and k'd have to be 
multiphed by 1/ \/2 or z/ \/2 for the cases that ui is either the scalar state or the pseudo-scalar 

^ In Eqs. (2.5) and (2.7) of [7,], the term q should correctly read uJq, as shown in Eq. (|A5|) here. 
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state, respectively. The total width is given by 

d^p, d^Pr S{Ei^ -E,-Er- Eo,) ,2 



16mt^(27r)5y Eb E^ 



Eo^ 



fi\ > 



with 






= 




= E 




i=l,2 



(B3) 



(B4) 



Tr- 



+ 2Re 



Tr 



12 



2 _ ^2 



(B5) 



+ 2m|.(ifc*,nzrir+i/*iri^rir)p.-p. 

- 4m,m-+Re(A;ri/ri*) (141' + iPh ' Pb) ■ Pb 

- ArrirmbRe (klik-il'^Jll^ ((p^^ - pb) ■ (pi, -pb)+ml+^ , 

Tri2 = 2 (/c*;^/c5';^/c2i/c2i + ^11^11^21^21) ^ 

(2(Pti - Pb) • Pr (Pti - Pb) • Pb - (Pti - Pb)' Pb • Pr) 

+ 2 m^+m^+ (^*nin^2i4*i + ^n^n 4*i^2i) Pr • Pb 

+ 2m6m-+ (Mi^n^2i4*i + ^ii^n ^21^2*1) (Pti -Pfe) "Pr 

+ 2m-+m6 (^ii^rr4*i^2i + ^n^n^2i4l) (Pf, - Pb) ■ Pr 

- 2m,m^+ (^*n^n ^21^2*1 + ^n^n^2i4*i) {Ph - Pb) ■ Pb 



(B6) 



- 2m^+m, [KiKmikt, + /ii/ri4*i^2ij (Pf^ -Po) -Pb 

- 2m,m6 (k\,kim,lt, + /*n/ri*^2i4*i) (Pt, -Pft) • (Pti -Pft) 

- 2m-+m,m6m-+(A;JiZrr4*i^2i + ^n^n^2i^2*i) ■ 



(B7) 
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